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Estimating the age of the Universe is an old problem. Rapid progress 
in observational cosmology in recent years has led to more accurate val- 
ues of the fundamental parameters. The current most popular model is 
a flat Universe, with about 30% of the critical density in the form of mat- 
ter (baryonic and non-baryonic) and 70% in the form of 'dark energy'. 
These two densities, together with the Hubble constant (estimated to be 
about 70 km s^^ Mpc^), constrain the age of the Universe to be approx- 
imately 13 Gyr. This expansion age is uncomfortably close to the age of 
the oldest globular clusters (approximately 12.5 Gyr), in particular if they 
formed relatively recently. We review here proposed models of globular 
cluster formation and point out possible conflicts with cosmology. 

In a Big Bang model with a cosmological constant A the age of the Uni- 
verse is determined by three present-epoch parameters: the Hubble con- 
stant Ho = 100 /i km s^^ Mpc^^ = (9.78 Gyr)^^ h, the mass density pa- 
rameter iljn/ and the scaled cosmological constant = 
cosmological constant tends to stretch the age of the Universe, but if the 
dark energy were in the form of 'quintessence' which decays with time, 
then tn would be reduced. For this reason (and to abbreviate our discus- 
sion) we focus on the 'optimistic' case of constant A. For a flat universe 
(Qjn -I- r^A = 1)/ supported by the recent Cosmic Microwave Background 
(CMB) experiments, the age depends only on two free parameters and is 
well approximated by: 

2 

This gives us an insight to the way the errors propagate in the determina- 
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tion of the cosmic age: 

-ii « — ^ + 0.3—^ . (2) 

This shows that the fractional error in Hq is three times as important as the 
fractional error in 0,^^- Typically the quoted error in Hq is 10%@, while the 
recent range of is of order 50%, so the expected fractional error in age is 
25% (e.g. for in ~ 13 Gyr, Atu^^ Gyr). 

Numerous studies O H O have compared and combined in a self -consistent 
way the most powerful cosmic probes: the CMB, galaxy redshift surveys, 
galaxy cluster number counts, type la supernovae, and galaxy peculiar ve- 
locities. These studies indicate that we live in a flat accelerating Universe, 
dominated by cold dark matter (CDM) and 'dark energy' (the cosmologi- 
cal constant A or some generalization such as 'quintessence' S). More pre- 
cisely, the data are consistent with a A-CDM model with Q.^ = 1 — ^0.3 
and h ^ 0.7, which corresponds to an expansion age tn = 13.5 Gyr. 

While the above A-CDM model is currently very popular, there is no 
simple theoretical explanation for the fortuitous near-equality of the present- 
epoch matter density Q,ni and 'dark energy' density JIa, nor for the true na- 
ture of these components. As a further diagnostic, we would like to revive 
an old conundrum: is the age of the Universe compatible with the ages of 
the oldest objects within it? For consistency, the ages of globular clusters 
{tec) must satisfy the following relation with the epoch of their formation 
(tf) and the age of the Universe: 

tH = tf + tGC- (3) 

The age of the oldest GC is estimated to be tac = 12.5 ± 1.2 Gyr0. Ra- 
dioactive dating of a very metal-poor star in the Galaxy (using ^^^U) gives 
a similar age of 12.5 ± 3 GyrS. Also, recently discovered Lya emitters IB 

and Ly break galaxies® at redshift z ~ 2 — 3 show an already evolved 
stellar population as old as 1 Gyr, comparable to the expansion age at that 
epoch. 

We see that the ages of old objects might be uncomfortably close to the 
age of expansion. It is commonly assumed that tf < 2 GyrW, and hence 
tf is neglected in eq. @. However, we point out here that is model 
dependent, and different models predict a wide range of ages. 

Formation of globular clusters in the Milky Way intimately relates to 
the formation scenario of the Galaxy itself. For decades, two apparently 
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conflicting models dominated the thought on Galaxy formation: 'mono- 
lithic collapse' versus highly fragmented star formation@. The real- 
ization that galaxies are embedded into dark matter halos0 has led to 
a constructive synthesis with cosmology. According to the CDM model, 
galaxies form as a result of gravitational growth and interaction of primor- 
dial fluctuations. Small objects collapse first and merge into larger systems, 
extending the hierarchy to progressively higher masses. The present halo 
of the Milky Way formed in dozens of mergers of smaller progenitors. 

Our understanding of star formation has also improved dramatically. 
Theoretical models still have many shortcomings 0, but the observational 
picture is becoming clearer. It looks that most stars form in clusters and 
associations of various sizes. The hierarchy of cluster masses ranges from 
the young OB associations to the old massive globular clusters, with no 
special scale between 10 and 10^ M©: dN/dM cx M"", a = 1.5 - 2IIIP1. 
On a global scale, the efficiency of globular cluster formation remains low. 

In large and small elliptical galaxies, McLaughlin finds the same ratio of 
the mass of globular clusters to the total baryonic (stellar + gaseous) mass 
of their host galaxy, ecc = Mqc/M^^, ^ 0.0026 ± 0.0005. 

The observed distribution within galaxies and cluster mass function 
dN/dM differ from the initial ones due to dynamical evolution. Small- 
mass clusters {M < 10^ M0) are gradually destroyed by stellar two-body 
relaxation and tidal interaction with the Galaxy0ll5 |2C]|^|^ while very 
massive clusters sink to the center via dynamical friction. However, most 
of the high-mass clusters are essentially unaffected by the evolution, and 
therefore, preserve the shape of the initial mass function. 

The suggested scenarios of globular cluster formation can be grouped 
into four main types, in decreasing order of the redshift of formation zj. 
For reference, in our assumed cosmology (Jim = 0.3, JIa = 0.7, h = 0.7), 
redshifts Zf = 7, 3, 1 correspond to the formation epochs tj = 0.8, 2.1 and 
5.8 Gyr, respectively. 

• Cosmological objects 

Peebles & Dicke@ were the first to propose that globular clusters are 
cosmological objects formed soon after recombination, with masses 
related to the Jeans mass (smallest mass of gas clouds able to col- 
lapse). However, including dark matter in the calculation would re- 
duce the mass of proto-globular clouds. Also, the expected efficiency 
of the conversion of gas into star clusters would be much higher than 
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the observed Egc- 

Hierarchical population (zj ~ 7 — 10) 

A hierarchical formation is more promising within small galaxies - 
progenitors of the Milky Way. In progenitors just massive enough 
to be cooled by atomic hydrogen, the gas collapses into a small disk 
where gravitational instability causes fragmentation into clouds and 
formation of star clusters of progressively larger sizes. Each progeni- 
tor galaxy may form a few globular clusters, still remaining very gas- 
rich before merging into the Galactic halo. Similar ideas have been 



put forward in the literatureE3'E3c3E3. See Figure || for more de- 
tail. 

Large galactic halo {zf 1 — 3) 

Globular cluster formation as a result of thermal instability in hot 
gas at virial temperature Tvir 10^ K has been suggested by Fall 
& ReesiS (variations of this idea include '^'^). Dense clouds of 



metal-poor gas would cool to Tc ~ 10^ K, confined by the pressure 
of diffuse gas, and would have the right Jeans mass to form massive 
globular clusters. The epoch of formation is constrained to be z/ < 3, 
when halos with large enough Tvir have been assembled. On the other 
hand, massive gaseous clouds must not have disrupted the thin disk 
of the Galaxy (7-10 Gyr old), which demands zj > 1. 

Mergers of disk galaxies {zf < 1) 

Finally, the observations of young massive star clusters in colliding 



galaxies prompted a "merger" scenario ffl. In this model, a merger 
of two spiral galaxies leads to a burst of star formation producing a 
large population of metal-rich clusters, in addition to the older popu- 
lation associated with the original spirals. We include this model for 
completeness, although it has been originally designed to account for 
the bimodal metallicity distribution of globular clusters in elliptical 
galaxies and predicts that most of the clusters form fairly recently. 

In the table we provide a comparison of the formation models, as they 
score against the observational properties of the Galactic globular clusters. 
The mark / shows the model meets the observational constraint, ? indicates 
possible but not certain agreement, and X denotes clear disagreement with 
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the data. Hierarchical model seems to fare better in this comparison, but 
scenarios with later formation are still possible. 

Figure ^ illustrates eq. @ for the three possible values of the formation 
redshift. Early formation (zj = 7) is favorite, zj = 3 is marginally consistent 
with the errors of Hq and globular cluster ages, but any lower redshift is 
ruled out in the flat Universe model with Qm = 0.3. 

Figure § shows the constraints on Hq and flm, assuming that the abso- 
lute ages of globular clusters are fixed. Again, the current 'best-fit' values 
are consistent only with Zf >3. 

We have thus compared cosmic ages as estimated by researchers in 
three different (somewhat unrelated) fields of astrophysics: the age of ex- 
pansion of the Universe, the age of globular clusters, and the epoch of their 
formation. We find that the three are consistent (within the errors) only 
over a relatively small parameter space, and hence there might be a prob- 
lem with any one of the estimates. This leads to several possible conclu- 
sions: 

(i) If the currently popular cosmological model (ilm = 0.3, Qa = 0.7,h = 
0.7) is correct and the age estimates of globular clusters are reliable, then 
one can put an extra constraint on the models of globular cluster formation, 
such that Zf >3. 

(ii) If we trust the cosmology and wish to have a late globular cluster 
formation, then a revision of the globular cluster ages (which are somewhat 
model-dependent) is required. 

(iii) If future research suggested that even old globular clusters formed 
at zj < 3 and the age estimates of globular clusters are valid, then the 
above comparison indicates a problem for the 'standard' A-CDM model. 
This is because the time-lag between the Big Bang and the era of globular 
cluster formation would then be > 2.1 Gyr, larger than the estimated errors. 
The constraint is even more severe if the "dark energy" is time-dependent 
quintessence rather than a constant A. Our overall conclusion is that there 
is a strong cosmological (as well as astrophysical) motivation for firming 
up our understanding of globular cluster formation. 

Correspondence and requests for materials should be addressed to O. Y.G. 
(e-mail: ognedin@ast.cam.ac.uk). 
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Table 1: SCOREBOARD OF GLOBULAR CLUSTER MODELS 



Observational properties 


P&D 


Hierarchy 


F&R 


Merger 


1. Narrow range of metallicity within 










individual dusters ((5[Fe/H] 0.1) 


/ 


/ 
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2. Average metallicity Z ~ 0.03 Zq 


X" 


/ 


?6 


X^ 


3. Spherical geometry of 










globular cluster system 


/ 


/ 


?-X'^ 


X^ 


4. Low efficiency ^qc 


X/ 


/ 


/ 


79 


5. Spread of relative ages of 










oldest clusters ~ 2 Gyr 


X 




/ 


/ 


6. Young massive clusters have 










a broad mass function 


X* 


/ 


?i 


V 


7. Globulars have similar properties 










regardless of the size or morphology 










of the host galaxy 


/ 


/ 


? 
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Overall*^ 


X 


/ 


? 


X 



a in this model globular clusters form immediately after 

recombination when the intergalactic gas has few metals; 

self-enrichment is unlikely as globular clusters must 

collapse quickly enough to remain gravitationally-bound 
b also likely to have lower metallicity; pre-enriched gas 

would have cooled to a temperature lower than 10^ K 
c remaining gas in normal spiral galaxies is already enriched, 

Z ~ 0.1 — 1 overproducing metallicity in new GCs 
d spherical geometry is assumed, but cloud fragmentation is 

likely to proceed in highest density regions, i.e. in the disk 
e gas in mergers is contained within the orbital plane, 

and new GCs must have a highly flattened distribution 
/ this model assumes almost 100% efficiency of gas conversion 
g observations suggest a possibly higher efficiency 
h age spread could be attributed to different time of virialization 

of the progenitors of different mass and the statistical variance 

of the progenitors of similar mass, creating a range of 

redshifts of formation 
i this model reproduces only the median mass scale, M ~ 10^ Mq 
j in these models cluster mass would depend on the virial temperature 

of the halo and the metallicity of cooling gas 
k note that these marks apply only to the formation of globular 

clusters in the Galaxy. 
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Figure 1: In a hierarchical scenario that we propose, the oldest globular 
clusters form in small progenitor galaxies with the mass of the first-ranked 
cluster proportional to the amount of gas supply in the progenitor, fbM. 
Massive clusters are almost unaffected by dynamical evolution and serve 
as indicators of the initial distribution. The plot compares the mass func- 
tion of large globular clusters (asterisks; data from pl|]), renormalized by a 
factor {fb ecc)"^ ~ 3000, with the cumulative number of virialized progen- 
itors of the present Milky Way halo. For three different redshifts, A^(> M) 
is calculated using the extended Press-Schechter formalism©. The two 
distributions match at z « 7, so that the first clusters might have formed 
soon thereafter. The virialization epoch of the progenitors is expected to 
have a statistical spread, which produces a corresponding spread of globu- 
lar cluster ages, in agreement with observations. If the efficiency of cluster 
formation eoc is not universal, i.e. significantly higher in some progenitors 
and lower in others, the distribution may shift towards higher redshift. For 
instance, filled circles demonstrate the case where 10% of the progenitors 
form clusters with the efficiency 10 ecc- The expected epoch of formation of 
oldest clusters thus follows to be z/ « 7 — 10. 
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Figure 2: The age of expansion vs the age of globular clusters, as a function 
of the formation redshift Zf (eq. ^. Shaded box contains the current values 
of tn and tec within the observational errors. In a flat Universe with = 
l-flm = 0.7, the age of expansion is tn = ^Hq^ J^a^^/^ ln[(l + riA^/2)(l - 
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Figure 3: The Hubble constant Hq vs matter density Qm, required to give 
the age of expansion consistent with the redshift of globular cluster forma- 
tion Zf, assuming globular cluster ages are fixed at tec = 12.5 Gyr and the 
Universe is flat (cf Fig. ^. Dashed line shows the limiting case tn = tac- 
Shaded box contains the 'best-fit' values. 



